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ABSTRACT 

All-sky data from the Planck survey and the Meta-Catalogue of X-ray detected Clusters of galaxies (MCXC) are combined to investigate the 
relationship between the thermal Sunyaev-Zeldovich (SZ) signal and X-ray luminosity. The sample comprises ~ 1600 X-ray clusters with redshifts 
up to ~ 1 and spans a wide range in X-ray luminosity. The SZ signal is extracted for each object individually, and the statistical significance of the 
measurement is maximised by averaging the SZ signal in bins of X-ray luminosity, total mass, or redshift. The SZ signal is detected at very high 
significance over more than two decades in X-ray luminosity (IO*'erg/s < Ljoo^fe)"''' S 2 x 10*'erg/s). The relation between intrinsic SZ signal 
and X-ray luminosity is investigated and the measured SZ signal is compared to values predicted from X-ray data. Planck measurements and 
X-ray based predictions are found to be in excellent agreement over the whole explored luminosity range. No significant deviation from standard 
evolution of the scaling relations is detected. For the first time the intrinsic scatter in the scaling relation between SZ signal and X-ray luminosity 
is measured and found to be consistent with the one in the luminosity - mass relation from X-ray studies. There is no evidence of any deficit in 
SZ signal strength in Planck data relative to expectations from the X-ray properties of clusters, underlining the robustness and consistency of our 
overall view of intra-cluster medium properties. 
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1. Introduction 

Clusters of galaxies are filled with a hot, ionised, intra-cluster 
medium (ICM) visible both in the X-ray band via thermal 
bremsstrahlung and from its distortion of the cosmic microwave 
background (CMB) from inverse Compton scattering, i.e., the 
Sunyaev-Zeldovich effect (Sunyaev & Zeldovich 1970, 1972, 
SZ, hereafter). The SZ signal can be divided into a kinetic SZ 
and a thermal SZ effect, originating from bulk and thermal mo- 
tions of ICM electrons, respectively. Since the kinetic SZ is 
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a second-order effect, we only consider the thermal SZ effect. 
Because of the different scaling of SZ and X-ray fluxes with 
electron density and temperature, SZ and X-ray observations 
are highly complementary. The combination of information from 
these two types of observations is a powerful one for cosmologi- 
cal studies, as well as for improving our understanding of cluster 
physics (see Birkinshaw 1999, for a review). In this framework, 
it is paramount to investigate to what degree the ICM properties 
inferred from SZ and X-ray data agree. 

Unfortunately, there is currently no consensus on whether 
predictions for the SZ signal based on ICM properties derived 
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from X-ray observation agree with direct SZ observations, ham- 
pering our understanding of the involved physics. Lieu et al. 
(2006) find evidence of a weaker SZ signal in the 3-year WMAP 
data than expected from ROSAT observations for 3 1 X-ray clus- 
ters. Bielby & Shanks (2007) reach similar conclusions using 
the same WMAP data and ROSAT sample and the additional 
Chandra data for 38 clusters. Conversely, Afshordi et al. (2007) 
find good agreement between the strength of the SZ signal in 
WMAP 3-year data and the X-ray properties of their sample of 
193 massive galaxy clusters. The last findings are supported fur- 
ther by the results of Atrio-Barandela et al. (2008), whose anal- 
ysis is based on the same SZ data and a larger sample of 661 
clusters. Diego & Partridge (2010) argue that a large contami- 
nation from point sources is needed to reconcile the SZ signal 
seen in the WMAP 5-year data with what is inferred from a large 
sample of ROSAT clusters. However, using the same SZ data and 
a slightly larger but similar sample of ROSAT clusters, Melin 
et al. (20 11) find good agreement between SZ signal and expec- 
tations. The latter finding is confirmed by the work by Andersson 
et al. (2011), where high quality Chandra data for 15 South Pole 
Telescope clusters are used. Finally, the WMAP 7-year data anal- 
ysis by Komatsu et al. (2011) argues for a deficit of SZ signal 
compared to expectations, especially at low masses. 

Improved understanding of this issue is clearly desired, since 
it would provide invaluable knowledge about clusters of galax- 
ies and aid in the interpretation and exploitation of SZ surveys 
such as the South Pole Telescope (SPT, Carlstrom et al. 2011) 
survey, the Atacama Cosmology Telescope (ACT, Fowler et al. 
2007), and Planck^ (Tauber et al. 2010). Since August 2009 the 
Planck satellite has been surveying the whole sky in nine fre- 
quency bands with high sensitivity and a relatively high spa- 
tial resolution. Planck data thus offer the unique opportunity to 
fully explore this heavily debated issue. As part of a series of 
papers on Planck early results on clusters of galaxies (Planck 
Collaboration 201 le,c,f,g,h), we present a study of the relation- 
ship between X-ray luminosity and SZ signal in the direction 
of ~ 1600 objects from the MCXC X-ray clusters compilation 
(Piffaretti et al. 2011) and demonstrate that there is excellent 
agreement between SZ signal and expectations from the X-ray 
properties of clusters. 

The paper is organised as follows. In Sect. 2 we briefly 
describe the Planck data used in the analysis and present the 
adopted X-ray sample. In Sect. 3 we present the baseline model 
used in the paper The model description is rather comprehen- 
sive because the model is also adopted in the companion pa- 
pers on Planck early results on clusters of galaxies (Planck 
Collaboration 201 le,c,f,h). Sect. 4 describes how the SZ signal 
is extracted from Planck frequency maps at the position of each 
MCXC cluster and how these are averaged in X-ray luminosity 
bins. Our results are presented in Sect. 5 and robustness tests are 
detailed in Sect. 6. Our findings are discussed and summarised 
in 7. 

When necessary we adopt a ACDM cosmology with //q - 
70 km/s/Mpc, Q.m = 0.3 and D.a - 0.7 throughout the paper 
The quantity E{z) is the ratio of the Hubble constant at redshift 
z to its present value. Ho, i.e., E{z)^ = f^m(l + z)^ + Ha. 



' Planck (http://www.esa.int/P/anc/:) is a project of the European 
Space Agency (ESA) with instruments provided by two scientific con- 
sortia funded by ESA member states (in particular the lead countries 
France and Italy), with contributions from NASA (USA) and telescope 
reflectors provided by a collaboration between ESA and a scientific con- 
sortium led and funded by Denmark. 



The total cluster mass M500 is defined as the mass within the 
radius R500 within which the mean mass density is 500 times 
the critical density of the universe, Pa-itiz), at the cluster red- 
shift: M500 = 5^Pcrit(z) 500/?5g(). We adopt an overdensity of 
500 since Rsqq encloses a substantial fraction of the total viri- 
alised mass of the system while being the largest radius probed 
in current X-ray observations of large samples of galaxy clusters. 

The SZ signal is characterised by I500 defined as 
^a'^^-' ~ (crj / ms^c') J PdV , where Da{z) is the angular dis- 
tance to a system at redshift z, crj is the Thomson cross-section, 
c the speed of light, nig the electron rest mass, P = n^kTe the 
pressure, defined as the product of the electron number density 
and temperature and the integration is performed over the sphere 
of radius R500. The quantity F500 is proportional to the apparent 
magnitude of the SZ signal and Y^oo is the spherically in- 
tegrated Compton parameter, which, for simplicity, will be re- 
ferred to as SZ signal or intrinsic SZ signal in the remainder of 
the paper All quoted X-ray luminosities are cluster rest frame 
luminosities, converted to the [0.1-2.4] keV band. 

2. Data 

In the following subsections we present the Planck data and the 
X-ray cluster sample used in our analysis. In order to avoid con- 
tamination from galactic sources in the Planck data we exclude 
the galactic plane: | b \< 14deg from the maps. In addition, we 
exclude clusters located less than 1 .5x beam full width half max- 
imum (FWHM) from point sources detected at more than 10 cr in 
any of the single frequency Planck maps, because such sources 
can strongly affect SZ measurements. 

2.1. SZ data set 

Planck (Tauber et al. 2010; Planck Collaboration 2011a) is the 
third generation space mission to measure the anisotropy of the 
CMB. It observes the sky in nine frequency bands covering 30- 
857 GHz with high sensitivity and angular resolution from 31' 
to 5'. The Low Frequency Instrument LFI; (Mandolesi et al. 
2010; Bersanelli et al. 2010; Mennella et al. 2011) covers the 
30, 44, and 70 GHz bands with amplifiers cooled to 20 K. The 
High Frequency Instrument (HFI; Lamarre et al. 2010; Planck 
HFI Core Team 2011a) covers the 100, 143, 217, 353, 545, and 
857 GHz bands with bolometers cooled to 0. 1 K. Polarisation is 
measured in all but the highest two bands (Leahy et al. 2010; 
Rosset et al. 2010). A combination of radiative cooling and three 
mechanical coolers produces the temperatures needed for the 
detectors and optics (Planck Collaboration 2011b). Two Data 
Processing Centers (DPCs) check and calibrate the data and 
make maps of the sky (Planck HFI Core Team 201 lb; Zacchei et 
al. 201 1). Planck's sensitivity, angular resolution, and frequency 
coverage make it a powerful instrument for galactic and extra- 
galactic astrophysics as well as cosmology. Early astrophysics 
results are given in Planck Collaboration, 201 le-x. 

In this paper, we use only the six temperature channel maps 
of HFI (100, 143, 217, 353, 545 and 857 GHz), corresponding 
to (slightly more than) the first sky survey of Planck. Details of 
how these maps are produced can be found in Planck HFI Core 
Team (201 lb); Planck Collaboration (201 le). At this early stage 
of the Planck SZ analysis adding the LFI channel maps does not 
bring significant improvements to our results. We use the full 
resolution maps at HEALPix- nside=2048 (pixel size 1 .72') and 
we assume that beams are adequately described by symmetric 

- http://healpix.jpl.nasa.gov 
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Table 1. Values of the beam full width half maximum assumed for 
each of the six channel maps of HFI. 



frequency [GHz] 


100 


143 


217 


353 


545 


857 


FWHM ['] 


9.53 


7.08 


4.71 


4.50 


4.72 


4.42 


FWHM error ['] 


0.10 


0.12 


0.17 


0.14 


0.21 


0.28 



Gaussians with FWHM as given in Table 1 . Uncertainties in our 
results due to beam coiTections, map calibrations and uncertain- 
ties in bandpasses are small, as shown in Sect. 6 below. 

2.2. X-ray data set 

The cluster sample adopted in our analysis, the MCXC (Meta- 
Catalogue of X-ray detected Clusters of galaxies), is presented 
in detail in Piffaretti et al. (2011). The information provided by 
all publicly available ROSAT All Sky Survey -based (NORAS: 
Bohringer et al. (2000), REFLEX: Bohringer et al. (2004), 
BCS: EbeHng et al. (1998, 2000), SGP: Cruddace et al. (2002), 
NEP: Henry et al. (2006), MACS: Ebeling et al. (2007, 2010), 
and CIZA: Ebeling et al. (2002); Kocevski et al. (2007)) and 
serendipitous (160SD: Mullis et al. (2003), 400SD: Burenin 
et al. (2007), SHARC: Romer et al. (2000); Burke et al. (2003), 
WARPS: Perlman et al. (2002); Horner et al. (2008), and EMSS: 
Gioia & Luppino (1994); Henry (2004)) cluster catalogues was 
systematically homogenised and duplicate entries were carefully 
handled, yielding a large catalogue of approximately 1800 clus- 
ters. For each cluster the MCXC provides, among other quanti- 
ties, coordinates, redshifts, and luminosities. The latter are cen- 
tral to the MCXC and to our analysis because luminosity is the 
only available mass proxy for such a large number of X-ray clus- 
ters. For this reason we will focus here on how the cluster rest 
frame luminosities provided by the MCXC are computed. Other 
quantities such as total mass and cluster size will be discussed in 
Sect. 3 below, because they are more model dependent. 

In addition to being converted to the cosmology adopted in 
this paper and to the [0. 1 -2.4] keV band (the typical X-ray survey 
energy band), luminosities are converted to that for an overden- 
sity of 500 (see below). This allows us to minimise the scatter 
originating from the fact that publicly available catalogues pro- 
vide luminosity measurements within different apertures. 

Because cluster catalogues generally provide luminosities 
measured within some aperture or luminosities extrapolated up 
to large radii (total luminosities), the luminosities L500 provided 
by the MCXC were computed by converting the total luminosi- 
ties to Lsoo using a constant factor or, when aperture luminosities 
are available, by performing an iterative computation based on 
the R EXCESS mean gas density profile and Lsoo - Mgoo relation. 
The REXCESS L500 - Msm calibration is discussed in Sect. 3 be- 
low. While the comparison presented in Section 5.3 of Piffaretti 
et al. (2011) indicates that the differences between these two 
methods do not introduce any systematic bias, it is clear that 
the iteratively computed L500 are the most accurate. The itera- 
tive computation was possible for the NORAS/REFLEX, BCS, 
SHARC, and NEP catalogues. As shown in Piffaretti et al. 
(201 1) the luminosities L500 depend very weakly on the assumed 
L500 - M500 relation. Nevertheless, when exploring the different 
L500 - ^500 relations detailed below, we consistently recompute 
L500 using the relevant Lsoo - Msoo relation. 

In addition, we supplement the MCXC sample with z > 0.6 
cluster data in order to enlarge the redshift leverage. These ad- 



ditional high redshift clusters are collected from the literature 
by utilizing the X-Rays Clusters Database BAX ^ and perform- 
ing a thorough search in the literature. For these objects we col- 
lect coordinates, redshift, and X-ray luminosity. The luminosity 
values given in the literature are converted to the [0.1-2.4] keV 
band and adopted cosmology as done in Piffaretti et al. (2011). 
Because the available luminosities are derived under fairly dif- 
ferent assumptions (e.g., aperture radius, extrapolation methods, 
etc.) we do not attempt to homogenise them to the fiducial lu- 
minosity L500 as done in Piffaretti et al. (201 1). In almost all the 
cases the adopted luminosity is however either the total luminos- 
ity (i.e. extrapolated to large radii) or the directly the luminosity 
Lgoo. Given the fact that the difference between these is close to 
10 % and that uncertainties affecting luminosity measurement of 
high redshift clusters are much larger, we treat all luminosities 
as fiducial luminosity L500. 

The MCXC and z > 0.6 supplementary clusters located 
around the galactic plane (| b \< 14deg) or near bright point 
sources (> lOcr, distance < 1.5 x FWHM) are excluded from 
the analysis. The resulting sample comprises 1603 clusters, with 
845 clusters being members of the NORAS/REFLEX sample. 
There is a total of 33 supplementary z > 0.6 clusters located in 
the sky region selected in our analysis. 

In Fig. 1 we show luminosity and mass as a function of red- 
shift for the whole sample with the NORAS/REFLEX and sup- 
plementary z > 0.6 clusters displayed with different colours. 
The figure shows the different clustering in the L-z plane of 
RASS (mostly NORAS/REFLEX) and serendipitously discov- 
ered clusters. The NORAS/REFLEX clusters are central to our 
study for many reasons. First, being the most luminous and nu- 
merous, they are expected to yield the bulk of the SZ signal from 
known clusters. Second, their distribution in the sky is uniform: 
NORAS and REFLEX cover the northern and southern sky, re- 
spectively, with the galactic plane excluded {\ b \< 20deg). 
Finally, the NORAS/REFLEX sample was also used in Melin 
et al. (201 1) in an analysis equivalent to the one presented in this 
work but based on WMAP-Syi data. For these reasons we use 
NORAS/REFLEX clusters as control sample in our analysis. 

For simplicity, in the remainder of the paper the whole com- 
pilation of MCXC plus supplementary z > 0.6 clusters will be 
referred to as MCXC. The [0.1-2.4]keV luminosities of the clus- 
ters in our sample range from 1.53 x lO"**^ to 2.91 x lO'*^ erg s"', 
with a median luminosity of 0.95 x 10"*"* erg s"', and redshifts 
range from 0.0031 to 1.45. Notice that while the adopted sample 
essentially comprises all known X-ray clusters in the sky region 
of interest, its selection function is unknown. The latter issue and 
how we evaluate its impact on our results is discussed in Sect. 
3.1 below. 



3. The cluster model 

Our cluster model is based on the REXCESS, a sample expressly 
designed to measure the structural and scaling properties of the 
local X-ray cluster population by means of an unbiased, repre- 
sentative sampling in luminosity (Bohringer et al. 2007). The 
calibration of scaling relations and the average structural param- 
eters of such an X-ray selected sample is ideal because is not 
morphologically biased. Furthermore, the gas properties of the 
REXCESS clusters can be traced by XMM-Newton up to large 
cluster-centric distances, allowing robust measurements at an 
overdensity of 500. 



^ http://bax.ast.obs-mip.fr/ 
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Table 2. Values for the parameters of the adopted Lx- M relation. 



0.01 




10^ 



Fig. 1. Observed [0.1-2.4] keV band luminosities (right vertical axis) 
and inferred masses (left vertical axis) as a function of redshift. Shown 
are the MCXC (the NORAS/REFLEX control subsample in shown in 
red) and the supplementary clusters (blue dots). 

Since X-ray luminosity is the only available mass proxy for 
our large cluster sample, the most fundamental ingredient of the 
cluster model is the scaling relation between [0.1-2.4] keV band 
luminosity and total cluster mass, which is detailed in Sect. 3.1. 
Given a cluster redshift z, mass M500 and hence cluster size ^500, 
the universal pressure profile of Arnaud et al. (2010) is then used 
to predict the electronic pressure profile. This allows us to pre- 
dict Y500, the SZ signal integrated in a sphere of radius R^oo 
as summarised in Sect. 3.2. It is important to notice that the es- 
timated cluster size R^qo and the universal pressure profile are 
also assumed when extracting the SZ signal from Planck data as 
detailed in Sect. 4 below. 

In the following we describe the assumptions at the basis 
of om fiducial model and provide the adopted scaling laws. In 
addition, we also discuss how these assumptions are varied in 
order to investigate the robustness of our results. 



3.1. L 



500 



• Mgoo relation 



For a given [0.1-2.4] keV band luminosity L500 the total mass 
M500 is estimated adopting the REXCESS L500 - A^soo relation 
(Pratt et al. 2009): 



E(z) 



-7/3 



L5Q0 



1 044 erg s- 



M500 



3 X lO'^Mp 



(1) 



Because this relation has been calibrated using the low scatter 
X-ray mass proxy Yx (Kravtsov et al. 2006), the parameters 
Clm and q-lm depend on whether the slope of the underlying 
A/500 - Yx relation is assumed to be equal to the standard (self- 
similar) value of umyx - 3/5 or it is allowed to be a free pa- 
rameter, yielding qtmyx - 0.561 (see Eqs. 2 and 3 in Arnaud 
et al. 2010). In the reminder of the paper these two cases will 
be referred to as standard and empirical, respectively. Om fidu- 
cial model adopts the empirical case, which reflects the observed 
mass dependence of the gas mass fraction in galaxy clusters. It 
is thus fully observationally motivated. 

In addition to these two variations of the L500 - M500 re- 
lation, we also consider the impact of Malmquist bias on our 
analysis. To this end we perform our analysis using the L - M 
calibrations derived from REXCESS luminosity data corrected 
or uncorrected for the Malmquist bias. In the reminder of the 



QTMYx 


L-M 


log Clm 




C loi> L-IobM 


0.561 


REXCESS 


0.274 


1.64 


0.183 


0.561 


Intrinsic 


0.193 


1.76 


0.199 


3/5 


REXCESS 


0.295 


1.50 


0.183 


3/5 


Intrinsic 


0.215 


1.61 


0.199 



paper these two cases will be referred to as the intrinsic and 
REXCESS L500 - M500 relations, respectively. Notice that the 
difference between the intrinsic and REXCESS L500 - M500 rela- 
tions is very small at high luminosity (Pratt et al. 2009). Ideally, 
one should use the intrinsic L500 - A/500 relation and compute, 
for each sample used to construct the MCXC compilation, the 
observed L500 - Mgoo relation according to each survey selection 
function. Unfortunately this would be possible only for a small 
fraction of MCXC clusters because the individual selection func- 
tions of the samples used to construct it are extremely complex 
and, in most of the cases, not known or not available. Therefore 
we simply consider the intrinsic L500 - M500 relation as an ex- 
treme and illustrative case, since it is equivalent to assuming that 
selection effects of our X-ray sample are totally negligible. On 
the other hand, in particular for the NORAS/REFLEX control 
sample and at high luminosities, the REXCESS L500 - Afgoo re- 
lation is expected to be quite close to the one that would be 
observed in our sample. For these reasons, our fiducial model 
adopts the REXCESS L500 - A/soo relation and the intrinsic case 
is used to test the robustness of our results. 

These different choices result in four different calibrations 
of the L500 - A/soo relation. The corresponding best fitting pa- 
rameters are summarised in Table 2 (see also Arnaud et al. 
2010). Values are given for the fiducial case where the observed 
REXCESS Lx - Yx and M - Yx ais assumed as well as for 
the cases where these two assumptions are varied: i.e. intrinsic 
(Malmquist bias corrected) Lx - Yx relation and standard slope 
of the M-Yx relation ajviYx- The table also lists the intrinsic dis- 
persion in each relation, which we use to investigate the effect of 
scatter in the assumed mass-observable relation in our analysis. 

For a given L500 - A/500 relation we estimate, for each clus- 
ter in our sample, the total mass A/goo from its luminosity L500. 
When the latter is computed iteratively (see Sect. 2.2), the same 
L500 - Af500 relation is adopted for consistency. Finally, the clus- 
ter size or characteristic radius R^m is computed from its defini- 



tion: M500 = Upcm(z)500R^ 



500- 



3.2. The SZ signal 

As shown in Arnaud et al. (2010), if standard evolution is as- 
sumed, the average physical pressure profile of clusters can be 
described by 



Pir) = P500 



Afsoo 



3 X 1014 Mo / (c50Qx)y(l + (csooxrY-^ 



(2) 



with X = r//?5oo and ap - l/aMYx - 5/3. In the standard case we 
have ap = 0, while in the empirical case ap = 0.12. Notice that 
the most precise empirical description also takes into account a 
weak radial dependence of the exponent ap of the form ap = 
0. 12 -H a'pix). Here we neglect the radially dependent term since, 
as shown by Arnaud et al. (2010), it introduces a fuUy negligible 
correction. 
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Table 3. Parameters describing the shape of the pressure profile. 









t'soo 


r 


a 


P 


All 


0.561 


8.403 


1.177 


0.3081 


1.0510 


5.4905 


CC 


0.561 


3.249 


1.128 


0.7736 


1.2223 


5.4905 


MD 


0.561 


3.202 


1.083 


0.3798 


1.4063 


5.4905 


All 


3/5 


8.130 


1.156 


0.3292 


1.0620 


5.4807 



The characteristic pressure P500 is defined as 

P500 = 1.65 X lQ-'Eiz}'l' ^T!.. r keVcm-l (3) 

The set of parameters [Po,C5oo,y, a,jS] in Eq. 2 are constrained 
by fitting the REXCESS data and depend on the assumed slope 
of the M - Yx relation. In Table 3 we list the adopted best fit- 
ting values, which, as detailed in Sect. 4 below, are also used to 
optimise the SZ signal detection. Values are first given for the 
fiducial case where the observed M - Yx relation (with slope 
ayiY^-0-56Y) and the average profile of all REXCESS clusters 
are adopted. The values for the average cool-core (CC) and mor- 
phologically disturbed (MD) REXCESS profiles, that we use to 
estimate the uncertainties originating from deviations from the 
average profile (see Sect. 6), and the average profile derived as- 
suming a standard slope of the M - Yx relation (aMYx - 3/5), 
are also listed in the table. 

Because of the large number of free parameters, there is a 
strong parameter degeneracy and therefore a comparison of in- 
dividual parameters in Table 3 is meaningless. The parameters 
for the standard case are also listed in the table. 

The model allows us to compute the physical pressure profile 
as a function of mass /W500 and z and thus to obtain the Fsoo 
- Msoo relation by integration of P{r) in Eq. 2 within a sphere of 
radius The relation can be written as 

Di(z) Y500 = 2.925 X 10-5/(1) 

or, equivalently, 

F500 = 1.383 X lO-V(l) 

/ M500 „, .on j Da(z) . 2 

X E(zY' arcmin , (5) 

\3 X 1014 Mo/ \500Mpc/ 

where /(I) - 0.6145 and 7(1) = 0.6552 are numerical factors 
arising from volume integrals of the pressure profile in the em- 
pirical and standard slope case, respectively (see Arnaud et al. 
2010, for details). Combining Eqs. 1 and 4 gives 

Di(z)Ysoo = 2.925x 10-5/(1) 

. £|:^(^52L_)|*s«./3„p,., ,« 

Clm yiO^^^ergs '/J 

where q;ly = ^lm x Q'mYv- In the fiducial case aiy — 0.92, im- 
plying that Ysoo oc L^^q for our model predictions. The clus- 
ter model allows us to predict the volume integrated Compton 
parameter Y^qq for each individual cluster in our large X- 
ray cluster sample from its [0.1-2.4] keV band luminosity Lsoo- 
These X-ray based prediction can be computed for different as- 
sumptions about the underlying X-ray scaling relations (stan- 
dard/empirical and intrinsic /REXCESS cases) and compared 



with the observed SZ signal, whose measurement is detailed in 
the next section. 

To reiterate, om fiducial case assumes: empirical slope of 
M-Yx relation and REXCESS L500 - M500 relation. If not other- 
wise stated, in the remainder of the paper results for the fiducial 
case are presented and results obtained by varying the assump- 
tions are going to be compared to it in Sect. 6. 

For simplicity the cluster size and SZ signal for the 
MCXC clusters in Planck Collaboration (201 le) are provided 
in the standard M^qq - Yx slope case. As we show in Planck 
Collaboration (201 le) the effects of this on X-ray size and both 
predicted and observed SZ quantities for clusters in the Early 
Sunyaev-Zeldovich (ESZ) catalog are fully negligible with re- 
spect to the overall uncertainties. 

4. Extraction of the SZ signal 

4. 1 . Individual measurements 

The SZ signal is extracted for each cluster individually by cut- 
ting from each of the six HFI frequency maps 10° x 10° patches 
(pixels 1.72 arcmin) centred at the cluster position. The result- 
ing set of six HFI frequency patches is then used to extract the 
cluster signal by means of multifrequency matched filters (MMF, 
hereafter). The main features of the multifrequency matched fil- 
ters are summarised in Melin et al. (201 1) and more details can 
be found in Herranz et al. (2002) and Melin et al. (2006). 

The MMF algorithm optimally filters and combines the 
patches to estimate the SZ signal. It relies on an estimate of the 
noise auto- and cross-power-spectra from the patches. Working 
with sky patches centred at cluster positions allows us to get 
the best estimates of the local noise properties. The MMF also 
makes assumptions about the spatial and spectral characteristics 
of the cluster signal and the instrument. Our cluster model is de- 
scribed below, for the instrumental response we assumed sym- 
metric Gaussian beams with FWHM given in Table 1 . 

We determine a single quantity for each cluster from the 
Planck data, the normalisation of an assumed profile. All the 
parameters determining the profile location, shape and size are 
fixed using X-ray data. We use the profile shape described in 
Sect. 3 with C500, a, /3 and y fixed to the values given in Table 
3 and integrate along the line-of-sight to obtain a template for 
the cluster SZ signal. The integral is performed by considering 
a cylindrical volume and a cluster extent of 5 x Rsm along the 
line-of-sight. The exact choice of the latter is not relevant. The 
normalisation of the profile is fitted using data within a circular 
aperture of radius 5 x /?5oo for each system in our X-ray cluster 
catalogue, centring the filter on the X-ray position and fixing the 
cluster size to flgoo = /?50{)/£'a(-z)- Notice that the dependence of 
cluster size on X-ray luminosity is weak (Rsqo oc L^qq from Eq. 
1), implying that MMF measurements are expected to be rel- 
atively insensitive to the details of the underlying L500 - M500 
relation. 

The MMF method yields statistical SZ measurement errors 
CT; on individual meaurements. The statistical error includes un- 
certainties due to the instrument (beam, noise) and to the astro- 
physical contaminants (primary CMB, Galaxy, point sources). 
Obviously, it does not take into account the uncertainties on our 
X-ray priors and instrumental properties which will be studied 
in Sect. 6 

The same extraction method is used in Planck Collaboration 
(201 Ih) where the optical-SZ scaling relation with MaxBCG 
clusters (Koester et al. 2007) are investigated. There are only two 
differences. First, in this paper we use the X-ray scaling L500 - 
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Fig. 2. Left: Intrinsic SZ signal from a sphere of radius R^qo as a function of the X-ray luminosity for all the clusters in the sample individually. 
Error bars indicate the pure measurement uncertainties based on MMF noise estimates (statistical uncertainties). Red diamonds show the bin 
averaged values with thick and thin error bars indicating the statistical (not visible) and bootstrap uncertainties, respectively. Right: Zoom onto 
the scale indicated by the horizontal dotted lines in the left-hand panel. Red symbols and error bars as in left-hand panel. Green triangles (shifted 
towards lower X-ray luminosity values by 20% with respect to diamonds for clarity) show the result of the same analysis when the signal is 
estimated at random positions instead of true cluster positions. The associated thick error bars indicate the statistical uncertainties. 



Msoo to adapt our filters to the sizes of our clusters while we 
use the optical A^2oo - A^500 relation of Johnston et al. (2007) 
and Rozo et al. (2009) in the other paper Second, the MaxBCG 
catalogue includes ~ 14,000 clusters so we do not build a set 
of patches for each cluster individually. Instead, we divide the 
sphere into 504 overlapping patches (10° x 10°, pixel=1.72 ar- 
cmin) as in Melin et al. (2011). We also extract SZ signal for 
each cluster individually but the clusters are no longer located at 
the centre of the patch. 

Under the assumption that the shape of the adopted pro- 
file template corresponds to the true SZ signal, our extraction 
method allows us to convert the signal in a cylinder of aperture 
radius 5 x /?5oo to Fsoo, the SZ signal in a sphere of radius Rsqq. 
By definition the conversion factor is a constant factor for every 
cluster but depends on the assumed profile. The effect of the un- 
certainties on the assumed profile are discussed in Sect. 6 below. 

The intrinsic SZ signal is computed by taking into account 
the angular distance dependence of the observed signal and is 
expressed as (Da(z)/500Mpc)- Fsoo. This signal has units of 
arcmin^ as for the observed quantity, but its value differs from 
the intrinsic signal in units of Mpc' by a constant, redshift- 
independent factor Making such a conversion allows us to di- 
rectly compare our measurements with the predictions derived 
from the model detailed in Sect. 3 (see in particular Eq. 4 and 5). 
When a specific scaling relation is investigated, the SZ signal is 
appropriately scaled according to the adopted scaling relations 
presented in Sect. 3 (e.g., see the left-hand panel of Fig. 2). 

The SZ signal for all the clusters in our sample is shown as a 
function of the [0.1-2.4] keV band X-ray luminosity in the left- 
hand panel of Fig. 2. Assuming standard evolution the intrinsic 
quantities (Da(z)/500Mpc)- Fsoo E{zT~^^ are plotted as a func- 
tion ofL5oo£'(z)"^^l Thefi gure shows that Planck detects the SZ 
signal at high significance for a large fraction of the clusters. 

4.2. Binned SZ signal 

As shown in the left-hand panel of Fig. 2 the SZ signal is not 
measured at high significance for all of the clusters. In particu- 
lar, low luminosity objects are barely detected individually. We 



therefore take advantage of the large size of our sample and av- 
erage SZ measurements in X-ray luminosity, mass, or redshift 
bins. The bin average of the intrinsic SZ signal is defined as the 
weighted mean of the signal in the bin (with inverse variance 
weight, cr7^, scaled to the appropriate redshift or mass depen- 
dence depending on the studied scaling relation) and the associ- 
ated statistical errors are computed accordingly. The binning de- 
pends on the adopted relation and will be detailed in each case. 

In the left-hand panel of Fig. 2 the binned signal is overlaid 
on the individual measurements. In this case the SZ signal is 
averaged in logarithmically spaced luminosity bins. We merged 
the lowest four luminosity bins into a single bin to obtain a sig- 
nificant result. The statistical uncertainties, which are depicted 
by the thick error bars, are not visible in the figure and clearly 
underestimate the uncertainty on the binned values. 

A better estimate of the uncertainties in the binned values 
comes from an ensemble of 10,000 bootstrap realisations of the 
entire X-ray cluster catalogue. Each realisation is constructed 
from the original data set by random sampling with replacement, 
where all quantities of a given cluster are replaced by those of 
another cluster. Each realisation is analysed in the same way as 
the original catalogue and the standard deviation of the average 
signal in each bin is adopted as total error. Bootstrap uncertain- 
ties, which take into account both sampling and statistical uncer- 
tainties, are shown by the thin error bars in the left-hand panel of 
Fig. 2. A visual inspection of the figure indicates that the SZ sig- 
nal is detected at high significance over a wide luminosity range. 
The lack of clear detection at LsooE(zy^'^ < 0.05 x 10'*'*erg/s is 
due to the combined effect of low signal and small number of 
objects. In the companion paper Planck Collaboration (201 Ih) 
we explore this low luminosity (mass) range in more depth. The 
results of these two complementary analysis are summarised and 
discussed in Sec. 7 (see Fig. 10 and related discussion). 

The difference between statistical and bootstrap errors are 
rendered in more detail in Fig. 3, where relative bootstrap uncer- 
tainties (dot-dashed line) are compared to in-bin relative statis- 
tical errors (solid line). The figure shows that for L^ooEizy^^^ < 
10''^erg/s statistical uncertainties are dominant. This implies that 
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Fig. 3. Bin averaged relative statistical errors (solid line) and relative 
bootstrap errors (dot-dashed line) are shown as a function of X-ray lu- 
minosity. The numbers given in the legend indicate the number of ob- 
jects in each luminosity bin. For comparison, the scaled unweighted 
standard deviation (dashed line) is also shown. 



intrinsic scatter, which is discussed in more detail in Sect.5.3, 
can only be measured at higher luminosity. 

Fig. 3 also shows the quantity (1/ Va^) x {(t^^^IY) (dashed 
line), which is computed from the unweighted raw scatter o-j-aw, 
the bin average Y, and the number of clusters in the bin N. The 
difference between the latter and the relative bootstrap uncer- 
tainties in the low luminosity bins is due to the range of relative 
errors on individual measurements in a given bin. 

As a robustness check, we have undertaken the analysis a 
second time using random cluster positions but keeping all the 
properties of our sample (sizes, profile shape). The result is 
shown by the green triangles in the right-hand panel of Fig. 2 
and, as expected, is consistent with no detection of the SZ sig- 
nal. This demonstrative null test clearly shows the efficiency of 
the MMF to pull out the SZ signal from our cluster sample. 
Additional robustness test are discussed in Sect. 6 below. 



5. Results 

5.1. TheDlYsoo- 



■ L500 and D\ Ysqq - M500 relations 



The main results of our analysis are summarised in Fig. 4. In the 
left-hand panel of the figure the individual and luminosity binned 
Planck SZ signal measured at the location of MCXC clusters are 
shown as a function of luminosity together with the luminosity 
averaged model predictions. The latter are computed by averag- 
ing the model prediction for individual clusters (see Sect. 3) with 
the same weights as for the measured signal. Notice that SZ sig- 
nal and X-ray luminosity are intrinsic quantities and are scaled 
assuming standard evolution. The figure shows the high signif- 
icance of the SZ signal detection and the excellent agreement 
between measurements and model predictions. The agreement 
between Planck measurements and X-ray based predictions is 
rendered in more detail in the right-hand panel of Fig. 4 where 
the f /fl«cA:-to-model ratio is plotted. Taking into account the to- 
tal errors given by the bootstrap uncertainties (thin bars in the 
figure), the agreement is excellent over a wide luminosity range. 
We model the observed Y^qq - L500 relation shown in the left- 



Table 5. Bin averages of the Z)^ Yjoo - 
panel of Fig. 4. 



relation shown in the left 



L 


500 


Nr. 




i'500 


AKsoo 


AKsoo 


range 


Obj. 






statistical 


total 


0.100 


- 0.222 


152 


0.162 


0.037 


0.006 


0.009 


0.222 


-0.331 


130 


0.272 


0.093 


0.009 


0.012 


0.331 


- 0.493 


144 


0.419 


0.169 


0.010 


0.012 


0.493 


- 0.734 


175 


0.615 


0.254 


0.012 


0.021 


0.734 


- 1.094 


190 


0.894 


0.401 


0.013 


0.020 


1.094 


- 1.630 


177 


1.319 


0.616 


0.016 


0.041 


1.630 


- 2.429 


149 


1.931 


0.879 


0.022 


0.057 


2.429 


- 3.620 


121 


2.997 


1.521 


0.026 


0.130 


3.620 


-5.393 


100 


4.138 


2.356 


0.038 


0.142 


5.393 


- 8.036 


51 


6.572 


3.456 


0.076 


0.171 


8.036 - 


11.973 


26 


9.196 


5.342 


0.126 


0.359 


11.973 


- 17.840 


9 


14.345 


7.369 


0.236 


1.758 



hand panel of Fig. 4 by adopting a power law of the form 



J'snn — Y^i 



500 



1044erg/s 



Da{z) ] 
SOOMpcj 



(7) 



and directly fitting the individual points shown in the figure 
rather than the binned data points. We use a non-linear least- 
squares fit built on a gradient-expansion algorithm (the IDL 
curvefit function). In the fitting procedure, only the statistical 
errors given by the MMF are taken into account. The derived 
uncertainties on the best fitting parameters are quoted in Table 
4 as statistical errors. In addition, uncertainties on the best fit- 
ting parameters are estimated through the bootstrap procedure 
described in Sec. 4. Each bootstrap catalogue fit leads to a set 
of parameters whose standard deviation is quoted as the uncer- 
tainty on the best fitting parameters. Values are given for three 
different choices of priors as given in Table 4, where the best fit- 
ting parameters are listed. The table also provides the prediction 
of our X-ray based model for comparison. 

Fixing the slope and the redshift dependence of the Fgoo - 
L500 relation, the best fitting amplitude is 0.451 x 10"^ arcmin-, 
in agreement with the model prediction 0.428 x 10""^ arcmin- 
at 1.8cr. When keeping the redshift dependence of the relation 
fixed but leaving the slope of the relation free, we find agree- 
ment between best fitting and predicted slopes at better than Icr, 
while the amplitudes remain in agreement at 1.3o-. For maxi- 
mum usefulness and in particular to facilitate precise compar- 
isons with our findings, we provide, in Table 5, the data points 
shown in the left-hand panel of Fig. 4. Values are given for 
the quantities Lgoo = L^ooEizy^^^ in units of lO'*^ erg/s and 
?500 = YsooEiz)-^'^ (Da(z)/500Mpc)^ in units of 10"^ arcmin^ 
Both total (i.e., bootstrap) and statistical errors on ?5oo are also 
listed. 

For completeness, we also investigate the Fgoo - Mgoo 
relation, where the masses Mgoo are computed from the L500 - 
M500 relation given in Eq. 1. Following the same procedure as 
for the Z)^ Fgoo - ^500 relation, we fit individual points of the 

Ysm - M500 plane with 



Y5QO - j^500,M 



3 X lO'^Mp 



Eizf"- 



Da(z) \ 
500Mpcj 



(8) 



The same cases as for the Ysoo - L500 relation are considered 
and the best fitting parameters are provided in Table 6 along 



7 



Planck Collaboration: Planck early results X. 



O Plonck 
^ Model 



0.1 



■•.'V-,;*'s.??^-'. 



1.0 



10.0 



W,_2..],,v [10" erg/s] 




W'-2..]kev [10" erg/s] 



Fig. 4. Left: Scaling relation between Planck SZ measurements and X-ray luminosity for ~ 1600 MCXC clusters. Individual measurements are 
shown by the black dots and the corresponding bin averaged values by the red diamonds. Thick bars give the statistical errors, while the thin bars 
are bootstrap uncertainties. The bin-averaged SZ cluster signal expected from the X-ray based model is shown by the blue stars. The combination 
of the adopted 7500 - Af soo and L500 - M500 relations (Eq. 6) is shown by the dashed blue line while the red dot-dashed line shows the best 
fitting power-law to the data (Eq. 7 and Table 4). Right: Ratio between data and model bin averaged values shown in the left panel. Error bars are 
as in the left panel. 

Table 4. Best fitting parameters for the observed Z)^ F500 - Lsoo relation given in Eq. 7. Values are given for three different choices of priors and 
as predicted from X-rays for comparison. Both total errors from bootstrap resampling and statistical errors are quoted. 





5^500,L [lO^^arcmin^] 








0.451 ± 0.003 Stat [±0.013 tot] 


1.087 (fixed) 


2/3 (fixed) 


Pla?ick + MCXC 


0.447 ± 0.006 Stat [±0.015 tot] 


1.095 ± 0.008 Stat [±0.025 tot] 


2/3 (fixed) 




0.476 ± 0.006 Stat [±0.025 tot] 


1.087 (fixed) 


-0.007 ±0.1 54 Stat [±0.518 tot] 


X-ray prediction 


0.428 


1.09 


2/3 



Table 6. Best fitting parameters for the observed Ysqq - M500 relation given in Eq. 8. Values are given for three different choices of priors and 
as predicted from X-rays for comparison. Both total errors from bootstrap resampling and statistical errors are quoted. 





Y500M [10"' arcmin-] 








0.896 ± 0.007 Stat [±0.027 tot] 


1.783 (fixed) 


2/3 (fixed) 


Planck + MCXC 


0.892 ± 0.008 Stat [±0.025 tot] 


1.796 ± 0.014 Stat [±0.042 tot] 


2/3 (fixed) 




0.945 ±0.012 Stat [±0.049 tot] 


1.783 (fixed) 


-0.007 ±0.1 54 Stat [±0.518 tot] 


X-ray prediction 


0.850 


1.783 


2/3 



with the model prediction. Concerning the agreement between 
best fitting parameters and model predictions, the conclusions 
drawn for the Y^qo - L500 relation obviously apply also for 
the Dl Fsoo - M500. 

5.2. Redshift evolution 

We also considered the case where the redshift evolution of the 
scaling relations is allowed to differ from the standard expecta- 
tion. Using the simplest model (Eq. 7 or equivalently Eq. 8) we 
attempt to constrain the power law index j6l (or equivalently ^Sm)- 
We find that the measured SZ signal is consistent with standard 
evolution (see Table 4) and our constraints on any evolution are 
weak. Fig. 5 shows the measured and predicted, redshift binned, 
SZ signal, the expected standard redshift evolution, and the best 
fitting model. The figure shows that, although measurements and 
predictions agree quite well, the best fitting model is constrained 



primarily by the low redshift measurements. Possible future im- 
provements are discussed below in Sect. 7. 

5.3. Scatter in fhe D\ Y500 - L500 relation 

As discussed in Sect. 4.2, we find a clear indication of intrinsic 
scatter in our measurements of the Y^m - Lgoo relation. In 
this section we quantify this scatter and discuss how our mea- 
surement compares with expectations based on the representa- 
tive REXCESS sample (Arnaud et al. 2010) and the findings re- 
ported in the companion paper discussing high quality observa- 
tions of local clusters (Planck Collaboration 201 Ig). 

The intrinsic scatter cri„,j is computed in luminosity bins 
as the quadratic difference between the raw scatter cr,-a„ (see 
Sect. 4.2) and the statistical scatter expected from the statis- 
tical uncertainties, i.e. cr^ , = cr?,^„ - ai,„,. The latter is esti- 

' intr law sia.r 

mated by averaging the statistical uncertainties in a given bin. 
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Fig. 5. Bin averaged SZ signal from a sphere of radius Rsqq (Ksoo) scaled 
by the expected mass and angular distance dependence as a function of 
redshift. The Planck data (red diamonds) and the SZ cluster signal ex- 
pected from the X-ray based model (blue stars) are shown together with 
the expected standard redshift evolution (dahed line). The best fitting 
model is shown by the dot-dashed line and the Icr confidence region is 
limited by the dotted lines. Here M500 is given in units of 3 x IO'^Mq. 



A^"' Yj cr^, where N is the number of clusters in the 
bin. For a given luminosity bin, the uncertainty Ao-jnti on the esti- 
mated intrinsic scatter are evaluated by (Acrinu-)^ - o"fnti-(2 N (N- 
I)r'2:(l+(cr2/tr2^.,.))2. 

We find that intrinsic scatter can be measured only for 
LsooE{zy^^^ > 10^"^erg/s, because the statistical uncertainties at 
lower luminosities are close to the value of the raw scatter (see 
also Sect. 4.2). In a given bin with average signal Y, the resulting 
fractional intrinsic scatter (Tina /Y is shown in Fig. 6 along with 
the fractional raw and statistical scatters. The estimated intrinsic 
scatter is close to 40 - 50% and in agreement with the expec- 
tations given in Arnaud et al. (2010) (o-|ogY5„o = 0.184 + 0.024, 
the range of these values is indicated by the coarse-hatched re- 
gion in the figure). Notice that the intrinsic scatter reported in 
Arnaud et al. (2010) is computed for the R EXCESS sample and 
evaluated adopting XMM-Newton luminosities and a predicted 
SZ signal for individual objects based on the same model as- 
sumed here but relying on the mass proxy Yx- Therefore, the 
intrinsic scatter quoted in Arnaud et al. (2010) reflects the in- 
trinsic scatter in the underlying L500 - M500 relation. In Planck 
Collaboration (201 Ig), where a sample of clusters detected at 
high signal to noise in the Planck survey (the ESZ sample, see 
Planck Collaboration 201 le) and with high quality X-ray data 
from XMM-Newton is used, the intrinsic scatter in the Ysoo 
- L500 relation is found to be criogY,^,, = 0.143 + 0.016. These 
values are shown in Fig. 6 by the fine-hatched region. In Planck 
Collaboration (201 Ig) it is found that cool core clusters are re- 
sponsible for the vast majority of the scatter around the relation. 
Because the sample used in this study is X-ray selected, we ex- 
pect it to contain a higher fraction of cool core systems than 
in the ESZ subsample studied in Planck Collaboration (201 Ig). 
This implies that the scatter in the Ysqq - isoo relation mea- 
sured in our sample is expected to be higher than the one found 
in Planck Collaboration (201 Ig), as observed. Given the segre- 
gation of cool core systems in the Ysoa - ^soo reported in 
Planck Collaboration (201 Ig), we investigate the link between 
the intrinsic scatter in the relation and cluster dynamical state 



Fig. 6. Fractional raw (dot-dashed blue line and triangles), statistical 
(dot-dot-dashed green line and plus signs), and intrinsic (dashed red 
line, diamonds, and error bars) scatter on the Z)^ Fjoo - L^oa relation. 
The coarse/fine-hatched regions corresponds to the 1 a uncertainties 
on the intrinsic scatter reported in Arnaud et al. (2010) and Planck 
Collaboration (201 Ig), respectively. 



using our large X-ray sample. To this end we compare Planck 
measurements and the X-ray based predictions (i.e., Eq. 6) for 
individual objects. In Fig. 7 we show the difference between 
Planck measurement and the X-ray based prediction in units of 
the measurement statistical error cr; (see Sect. 4) as a function of 
X-ray luminosity and investigate the largest outliers the figure, 
i.e. the most statistically significant outliers in the 1500 - ^500 
relation. 

Given the size of our sample we discuss individually only 
clusters for which measurement and prediction differ by more 
than 5 0-;. These are identified and information on their dynam- 
ical state is searched for in the literature. Information is based 
on the classification of Hudson et al. (2010) if not stated other- 
wise. We find 15 clusters with a predicted signal smaller than the 
Planck measurement by more than 5 cr;. Of these, six are known 
merging clusters: Coma, A2218 (Govoni et al. 2004), 1ES0657, 
A754 (Govoni et al. 2004), A2163 (Bourdin et al. 201 1), A0697 
(Girardi et al. 2006), six are classified as non-cool core clusters 
and may therefore be unrelaxed: A2219 (Allen & Fabian 1998), 
A2256, A2255, A0209 (Zhang et al. 2008), A2813 (Zhang et al. 
2008), A3404 (Pratt et al. 2009), and A3266 is a weak cool core 
cluster No information is available for the remaining clusters: 
A1132and A3 186. Conversely, there are 11 over-predicted clus- 
ters at 5 CT;. Of these five are strong cool core clusters: 2A0335, 
Zwl021.0-f0426 (Morandi et al. 2007), A3112, HerA (Bauer 
et al. 2005) , and A0780. No information is available for the 
remaining clusters: A689, ACOSllll, A3392, J1253.6-3931, 
J1958.2-301 1, and RXCJ0643.4-I-4214. Notice that the luminos- 
ity of A689 is likely to be overestimated by a large factor be- 
cause of point source contamination (Maughan 2008). In addi- 
tion, for A3 186 and ACOSllll the model prediction rely on 
the EMSS luminosity measurements given in Gioia & Luppino 
(1994), which might be unreliable. 

The high fraction of dynamically perturbed / cool-core clus- 
ters with largely under/over predicted SZ signal is confirmed 
when additional outliers at smaller cr; are searched. We find 65 
clusters with (isoo Planck - Kgoo Model) > 3 cTi and that for 46 
percent of them dynamical state information is available in the 
literature. Of the latter 26 (87 percent) are either known mergers 
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Fig. 7. Difference between the Planck measurement and the X-ray 
based prediction in units of the measurement statistical error ct; (pure 
measurement uncertainties based on MMF noise estimates) as a func- 
tion of X-ray luminosity. Labelled black points denote objects with a 
difference larger than 5 cr\ and are further discussed in the text. Clusters 
with SZ signal possibly contaminated by radio sources (see discussion 
in Sect. 6) are shown in red and labelled by their name. 

or non-cool core clusters, two are weak cool core clusters, and 
only one is classified as cool core cluster. We find 53 clusters 
with (F500 Planck- F500 Model) < -3 (Ti. For 45 percent of these 
clusters we are able to find information on their dynamical state 
in the literature and find that 96 percent of them are cool core 
clusters. These findings clearly suggest that the intrinsic scatter 
in the D\ Y^qq - L500 relation is linked to the cluster dynamical 
state, as also found in Planck Collaboration (201 Ig). 

6. Robustness of the results 

As in the other four Planck SZ papers (Planck Collaboration 
201 le,c,g,h), we test the robustness of our results for the effect 
of several instrumental, modelling, and astrophysical uncertain- 
ties. Tests common to all Planck SZ papers are discussed in de- 
tail in Sec. 6 of Planck Collaboration (201 le). Of these, calibra- 
tion and colour correction effects are relevant for our analysis. 
Calibration uncertainties are shown to propagate into very small 
uncertainties on SZ signal measurements (~ 2%) and colour cor- 
rection is found to be a ~ 3% effect for Planck bands. 

In the following we report on robustness tests aimed at com- 
pleting this investigation. We show that our results are robust 
with respect to the instrumental uncertainties, that they are in- 
sensitive to the finest details of our cluster modelling, and that 
they are unaffected by radio source contamination. We show that 
restricting the analysis to the reference homogeneous subsample 
of NORAS/REFLEX clusters leads to measurements fully com- 
patible with those we obtain for the whole sample. 

6. 1 . Beam effects 

The beam effects studied in Planck Collaboration (201 le) are 
further scrutinised by directly estimating their impact on our 
results. To this end, the whole analysis is redone by assuming 
different beam FWHM. For simplicity, we systematically in- 
crease/decrease the adopted beam FWHM for all channels si- 



multaneously by adding/removing the conservative uncertainties 
given in Table 1 from the fiducial beam FWHM values. We find 
that the binned SZ signal varies by at most 2% from the value 
computed using the fiducial beams FWHM. 

6.2. Modelling 

The effects of changes in the underlying X-ray based model 
on our results are investigated by repeating the full analysis as 
for the fiducial case, but by assuming the standard slope of the 
M^oo - Yx relation and/or the intrinsic Lsoo - Msoo relation (see 
Sect. 3). For simplicity, in the following we discuss results ob- 
tained by varying only one assumption at a time. We find that 
the effect resulting by varying both assumptions is equivalent to 
the sum of the effects obtained by varying the two assumptions 
separately. 

As expected from the weak dependence of cluster size R^oo 
on luminosity, the measured SZ signal is barely affected by these 
changes. If the standard slope case is adopted instead of the em- 
pirical one, the bin averaged SZ signal changes by less than a 
few percent at all luminosities and the same is found when the 
intrinsic L500 - M500 relation is adopted. The model predictions 
are of course more affected by changes in the assumed scaling 
relations. In Fig. 8 we contrast the Planck -to-model ratio ob- 
tained for the different cases. The figure shows that the assump- 
tion on the slope of the M500 - Yx relation has a fully negligible 
impact. On the other hand it shows that the intrinsic L500 - Mgoo 
relation is not compatible with our measurements (> 5cr dis- 
crepancy). This finding is not surprising given the fact that when 
adopting the intrinsic L500 - M500 relation one assumes that se- 
lection effects of our X-ray sample are negligible. Notice that 
the WMAP-Syi data used in the similar analysis by Melin et al. 
(201 1) did not have sufficient depth to come to this conclusion. 
Furthermore, the agreement of our results for the R EXCESS and 
intrinsic L500 - M500 relations at high luminosity confirms that 
Malmquist bias is small for very luminous objects. 

6.3. Intrinsic dispersion in the L500 -M500 relation 

The intrinsic dispersion in the L500 - M500 relation dominates the 
uncertainty on the clusters' size Rsqq in our analysis. We investi- 
gate how this propagates into the uncertainties on the binned SZ 
signal by means of a Monte Carlo (MC) analysis of 100 realisa- 
tions. We use the dispersion given in Table 2 and, for each reali- 
sation, we draw a random mass log M500 for each cluster from a 
Gaussian distribution with mean given by the Lgoo - M500 rela- 
tion and standard deviation criogL-iogM/o'M- For each realisation, 
we extract the signal with the new values of M500 (thus TJgoo). 
The standard deviation of the SZ signal for the 100 MC realisa- 
tions in a given luminosity bin is found to be at most ~ 3% of the 
signal. Hence, given the size of the total errors on the binned SZ 
signal (see Fig. 3) our conclusions are fully unaffected by this 
effect. 

6.4. Pressure profile 

Furthermore, we investigate how the uncertainties on the as- 
sumed pressure profile propagates into the uncertainties on the 
binned SZ signal. For simplicity, we only quantify the largest 
possible effect by redoing the analysis but adopting the pressure 
profile parameters for the cool-core and morphologically dis- 
turbed subsamples given in Table 3, i.e. we assume that all clus- 
ters in the sample are cool-core or morphologically disturbed. 
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Both of the two resulting sets of binned SZ signal deviate from 
the one derived assuming the universal pressure profile by ap- 
proximatively 8% in the lowest luminosity bin and decrease lin- 
early with increasing logLsoo, becoming approximatively 1 % in 
the highest luminosity bin. Furthermore the normalisation of Eq. 
4 changes by less than 3% if the average pressure profiles param- 
eters of cool-core and morphologically disturbed clusters given 
in Table 3 are adopted instead of the ones for the average profile, 
implying that our SZ signal predictions are robust. Considering 
the total errors and their trend with luminosity, we conclude that 
our findings are fully unaffected by the exact shape of the SZ 
template. 



6.5. X-ray sample 

Because of the reasons detailed in Sect. 2.2, we also repeated our 
analysis by considering the NORAS/REFLEX control sample 
and find results fully consistent with those derived for the full 
sample. The results are shown in terms of Planck -to-model ratio 
in Fig. 9. Notice that in this case the luminosity binning is chosen 
so as to be comparable with that in the WMAP-Syi analysis of 
Melin et al. (2011). The comparison between WMAP-Syi and 
Planck results is discussed in Sec. 7 below. 



6.6. Radio contamination 

In addition we investigated the effect of contamination by radio 
sources on our results. Most radio sources are expected to have a 
steep spectrum and hence they should not have significant fluxes 
at Planck frequencies. However, some sources will show up in 
Planck LFI and HFI channels if their radio flux is sufficiently 
high and/or their spectral index is near zero or positive. Extreme 
examples are the Virgo and Perseus clusters that host in their in- 
terior two of the brightest radio sources in the sky. In the ESZ 
sample (Planck Collaboration 20 lie) there are also a few exam- 
ples of clusters with moderate radio sources in their vicinity ( 1 
Jy or less in NVSS) and still significant signal at LFI (and even 
HFI) frequencies. To check for possible contamination we com- 
bine data from SUMSS (a catalog of radio sources at 0.85 GHz, 
Bock et al. 1999) and NVSS (a catalog of radio sources at 1.4 
GHz, Condon et al. 1998). We have looked at the positions of 
the clusters in our sample and searched for radio sources in a 
radius of 5 arcmin from the cluster centre. We find that 74 clus- 
ters have a radio source within this search radius in NVSS or 
SUMSS with a flux above 1 Jy. Among these, eight have fluxes 
larger than 10 Jy, two sources larger than 100 Jy and one is an 
extreme radio source with a flux larger than 1 KJy. 

As a robustness test, we investigate the impact of contamina- 
tion by radio sources on our results by excluding clusters hosting 
radio sources with fluxes larger than 1 or 5 Jy and comparing 
the results to those obtained for the full sample. Interestingly 
we find that, as expected, the individual SZ signal is on aver- 
age lower than the X-ray based predictions in clusters that are 
likely to be highly contaminated. This is shown in Fig. 7 where 
clusters associated with radio sources with fluxes larger than 5 
Jy are shown by the red symbols. However, given the very low 
fraction of possibly contaminated clusters, we find that bin aver- 
aged signal is fully unaffected when these are excluded from the 
analysis. 
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Fig. 8. Ratio of binned Planck data points to model for different model 
assumptions. The fiducial model (black diamonds) is shown together 
with results obtained by varying the underlying L500 - M^m relation 
from REXCESS to intrinsic (green plus signs), and by varying the 
slope of the underlying M500 - i'x relation from empirical to standard 
(red triangles). Thick bars give the statistical errors, while the thin bars 
are bootstrap uncertainties. 



7. Discussion and conclusions 

As part of a series of papers on Planck early results on clus- 
ters of galaxies (Planck Collaboration 201 le,c,f,g,h), we mea- 
sured the SZ signal in the direction of ~ 1600 objects from the 
MCXC (Meta-Catalogue of X-ray detected Clusters of galaxies, 
Piffaretti et al. 201 1, see Sect. 2.2) in Planck whole sky data (see 
Sect. 2. 1) and studied the relationship between X-ray luminosity 
and SZ signal strength. 

For each X-ray cluster in the sample the ampUtude of the 
SZ signal is fitted by fixing the cluster position and size to the 
X-ray values and assuming a template derived from the univer- 
sal pressure profile of Arnaud et al. (2010). The universal pres- 
sure profile was derived from high quality data from REXCESS. 
Recently, Sun et al. (201 1) found that the universal pressure pro- 
file also yields an excellent description of systems with lower 
luminosities than those probed with REXCESS. This implies 
that the adopted SZ template is suitable for the entire luminosity 
range explored in our work. 

The intrinsic SZ signal F500 is averaged in X-ray lumi- 
nosity bins to maximise the statistical significance. The signal 
is detected at high significance over the X-ray luminosity range 
10''^erg/s < LsooEizy'"^ < lO'^^erg/s (see Fig. 2). 

We find excellent agreement between observations and pre- 
dictions based on X-ray data, as shown in Fig. 4. Our results 
do not agree with the claim, based on a recent WMAP-lyr data 
analysis, that X-ray data over-predict the SZ signal (Komatsu 
et al. 2011). Due to the large size and homogeneous nature of 
the MCXC, and the exceptional internal consistency of our clus- 
ter model, we believe that our results are very robust. Moreover, 
as reported in Sect. 6, we show that our findings are insensitive 
to the details of our cluster modelling. Furthermore, we have 
shown that our results are robust against instrumental (calibra- 
tion, colour correction, beam) and astrophysical (radio contami- 
nation) uncertainties. 

Our results confirm to a higher significance the results of the 
analysis by Melin et al. (2011) based on WMAP-Syr data. This 
is shown in Fig. 9 where the data-to-model ratio as a function of 
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Fig. 9. Data-to-model ratio for Planck results for the full sample (black 
diamonds) and the NORAS/REFLEX control sample (green plus signs). 
The WMAP-Syi results for the NORAS/REFLEX by Melin et al. (201 1) 
are shown by the red triangles. Error bars are as in Fig. 4. 



luminosity is presented. Luminosity bins are chosen so as to be 
comparable to those of Melin et al. (201 1), and the Planck results 
are presented for the whole sample used in this work and also for 
the NORAS/REFLEX sample adopted in Melin et al. (2011). In 
addition to the good agreement between results from the two 
data sets, the figure shows that in the WMAP-Syi study by Melin 
et al. (2011) statistical uncertainties are dominant. As shown in 
Sect. 4.2, Planck data allows us to overcome this limitation and 
to investigate the intrinsic scatter in the scaling relation between 
intrinsic SZ signal Y^qo and X-ray luminosity L500 (see Sect. 
5.3). We find a ~ 40% intrinsic scatter in the F500 - ^500 
relation and show that it is linked to cluster dynamical state. 

The agreement between luminosity binned X-ray predictions 
and Planck measurements is reflected in the excellent accord 
between predicted scaling relation and best fitting power law 
model to the Fsoo - isoo relation. The power law fit, which 
is performed on individual data points, is compared by Planck 
Collaboration (201 Ig) to the calibration derived from a sample 
of galaxy clusters detected at high signal to noise in the Planck 
survey (the ESZ sample, see Planck Collaboration 201 le) and 
with high quality X-ray data from XMM-Newton. As discussed 
in Planck Collaboration (201 Ig) the slight differences between 
the two best fitting relations reflect the difference between the 
selection of the adopted samples. Indeed, the X-ray sample used 
in the present work is X-ray selected and therefore biased to- 
wards the cool core systems, while the sample used in Planck 
Collaboration (201 Ig) is SZ selected. 

As mentioned in Sect. 4.2, the luminosity range where we 
are not able to detect the SZ signal because of the small num- 
ber of low mass objects (see Fig. 2), is explored in Planck 
Collaboration (201 Ih). In the latter analysis we use the optical 
catalogue of ~ 14,000 MaxBCG clusters (Koester et al. 2007) 
and, in a fully similar way as done in this work, extract the op- 
tical richness binned SZ signal from Planck data. By combining 
these results with the X-ray luminosity of the MaxBCG clus- 
ters measured by Rykoff et al. (2008) by stacking RASS data, 
in Planck Collaboration (201 Ih) we derive the Ysm - ^500 
relation for the MaxBCG sample. This result is shown together 
with the one derived in the present paper in Fig. 10. The X-ray 
luminosity histograms shown in the top panel of the figure high- 
light the complementarity of the two analyses. The bottom panel 
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Fig. 10. Bottom panel: Comparison between our results (red di- 
amonds, as in left-hand panel Fig. 4) and those obtained by 
Planck Collaboration (201 Ih) (green triangles), where MaxBCG 
clusters are investigated. X-ray luminosities and associated er- 
ror bars for the MaxBCG clusters are based on the analysis of 
Rykoff et al. (2008). Vertical error bars are as in Fig. 4 and the 
X-ray prediction (i.e., Eq. 6) is shown by the dashed blue line. 
Top panel: X-ray luminosity histograms of the MCXC (red) and 
MaxBCG (green) samples. For the MCXC the width of the bars 
is equal to the luminosity bin width, while for the MaxBCG we 
adopt the horizontal error bar shown in the bottom panel. 

of the figure shows agreement between the results from the two 
data sets and, very importantly, that observations and predictions 
based on X-ray data agree over a very wide range in X-ray lumi- 
nosity. 

We investigate the evolution of the scaling relation and find 
it to be consistent with standard evolution. Although redshift 
binned measurements and predictions agree quite well over a 
wide redshift range (see Fig. 5), our constraints are weak be- 
cause the inferred best fitting model is almost completely con- 
strained by only the low redshift measurements. Given the rele- 
vance of SZ-selected samples for cosmological studies and the 
need for complementary X-ray observations for such studies (see 
Planck Collaboration 2011c, and discussion therein), improved 
understanding of the evolution of SZ-X-ray scaling relations is 
clearly desired. High quality data similar to those used in Planck 
Collaboration (201 Ig), but for higher redshift clusters will pro- 
vide tight constrains on evolution, in particular when newly SZ 
discovered clusters (see Planck Collaboration 201 le, and refer- 
ences therein) with high quality X-ray and optical data are in- 
cluded. 
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